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1 Introduction 



Three dimensional massive gravity model introduced by Bergshoeff, Hohm and Townsend (BHT) 
leads to some renewed interests in three dimensional higher curvature gravity, since this so-called 
new massive gravity (NMG) is composed of the standard Einstein-Hilbert term with a specific 
combination of scalar curvature square term and Ricci tensor square one [I]~[6]. One of the original 
motivation of this theory is the non-linear completion of Fierz-Pauli massive graviton theory or the 
realization of the massive gravity on three dimensional space. Though NMG has different nature, 
for instance parity even rather than odd, from topologically massive gravity (TMG) [^[8], these 
theories share common features in some aspects. One of the most interesting common aspects 
in these theories is the existence of (warped) AdS black hole solutions p][9]~p!2]. Though the 
incompatibility of the unitarity of massive graviton modes and AdS black hole solutions is not 
resolved, it was suggested that the special case may exist as a consistent quantum gravity fT3]. 

Aside from massive gravitons, AdS black hole solutions give us another motivation to study NMG 
in the viewpoint of the AdS /C FT correspondence, as was done extensively in TMG |21j^|30]. One 
direction along the AdS /C FT correspondence, a la holographic c-theorem, is the extension of 
NMG to even higher curvature gravities, named as extended NMG [31j[32j. In these cases, various 
AdS black hole solutions are found and the central charges of their conjectural dual CFT are 
identified [3l]~01]. On the contrary, there are different aspects in NMG from TMG, one of which 
is the existence of the new type black holes discovered by BHT in NMG which we will name new 
type black holes in the following These black holes exist only for a specific combination of 
parameters in the NMG Lagrangian. 

New type black holes are asymptotically AdS and contain two parameters even in the static 
case. Prom the standpoint of the usual Einstein-Hilbert gravity, the existence of an additional 
parameter is mysterious for static and spherically symmetric black holes (or axi-symmetric ones in 
the three dimensional case). Since NMG is a higher derivative theory, the existence of an additional 
parameter may be regarded as natural when one considers an initial value problem. Though the 
masses of these black holes are obtained as a conserved charges in the form of the combination of 
black hole parameters, and the validity of the formula is supported by the first law of black hole 
thermodynamics and also by the dynamical new type black hole solutions, it is still unknown what 
is the physical interpretation of additional parameter. Suppose that the additional parameter or 
another combination of two parameters different from mass combination is identified as a conserved 
charge, it is puzzling that the first law of black hole thermodynamics is already satisfied without a 
new conserved charge. Moreover, mass seems to be the unique conserved quantity for static new 
type black holes from the symmetry consideration. 

If the additional parameter is not a conserved charge, it is desirable to study non-conserved 
quantities relevant to black hole physics. One of these non-conserved dissipating quantities in black 
hole physics is quasi- normal mode(QNM) [45] ~ [47] (See references there in). Some information 
about black holes can be obtained by perturbing them and studying ensuing behaviors. Since black 
holes are not closed systems, there are dissipating modes when black holes are perturbed. The 
dissipating modes in black holes may be defined by modes such as only ingoing part at the horizon 
and only radiating one at the infinity. These modes with a specific boundary condition are called 
QNMs. The same term is used for modes of dissipating fields other than metric perturbation in 
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the black hole background. In this paper we study QNM frequencies, which are complex numbers, 
of scalar fields on the new type black holes. 



There are several suggestions how to define QNMs on black holes. One of them is to take the 
boundary conditions such as ingoing modes only at the horizon and outgoing modes only at the 
spatial infinity with complex frequencies. In the asymptotic AdS case, which is relevant to us, the 
Dirichlet boundary condition instead of outgoing one may be taken at the spatial infinity because 
of the existence of blowing up potential [H]. QNMs mean the decay of modes for the time elapse. 
Because of the complex nature of QNM frequencies, ingoing/outgoing boundary conditions lead to 
exponentially blowing behaviors of modes at the horizon and the spatial infinity for fixed time. By 
separation of variables, the radial part of field equations is given in the form of 



= 0. 



In the viewpoint of the radial equations, ingoing/outgoing boundary conditions mean taking the 
exponentially blowing part while ignoring exponentially vanishing one. This process may be am- 
biguous since we should take care of exponentially small parts which are usually difficult to control 
in approximate solutions or asymptotic series solutions of differential equations, though one may 
not confront such difficulties in exact solutions. 



To overcome the ambiguity in boundary conditions taken in the domain of real values, it was 
proposed to perform the analytic continuation of radial coordinate to the complex domain and 
to impose boundary conditions along Stokes lines [50] • Stokes lines are defined by the contour 
such that would-be vanishing or blowing terms are equally contributingj^] In our notation these 
lines can be defined by 3m{ujx) = 0. On these lines the boundary conditions can be imposed 
unambiguously. It may be useful to note that one needs some additional information about uj 
actually to draw Stokes lines on the complex plane of the radial coordinate. Along with these 
boundary conditions, a method is developed to obtain QNM frequencies using Stokes lines. This 
method is basically monodromy approach to QNMs, while one needs some modifications in the 
asymptotic AdS case [l9]~[5T]. 

In the asymptotic AdS case Stokes lines are not closed and so open lines are used to match 
approximate solutions. This Stokes line method gives us approximate results since it uses matchings 
among the approximate solutions through the interpolation. However, some important outcome 
in QNM frequencies can be obtained. For instance, the asymptotic QNM frequencies can be 
reliably identified, which is important in the area spectrum and entropy quantization through 
Hod's conjecture |52j|53j. In our case of new type black holes, one can see some clues, through 
QNM frequencies, about which combination of parameters different from mass is appropriate. 



This paper is organized as follows. In section 2, we give a brief review on the properties of new 
type black holes in NMG to fix our conventions. In section 3, we solve the scalar field equations 
on new type black hole backgrounds and obtain exact QNMs in some special cases, which include 
results in the case of the Bahados-Teitelboim-Zanelli (BTZ) black holes [5l]. In this section we have 
taken the simple form of boundary conditions for QNMs. In section 4, we obtain highly damped 
QNMs in some general case using Stokes line method. We conclude with summary and several 

^In mathematical literatures these contours are called anti-Stokes lines. 
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comments in the final section. In the appendixes, we present some mathematical formulae and 
numerical computation of QNM frequencies for some parameter domains. 



2 New Type Black Holes in NMG 



2.1 New massive gravity : brief review 

After Bergshoeff, Hohm and Townsend have introduced a higher curvature theory (NMG) as the 
non- linear completion of Fierz-Pauli massive graviton theory in three dimensions [I]~[l], of which 
Lagrangian consists of the standard Einstein-Hilbert term and a specific combination of scalar cur- 
vature square term and Ricci tensor square one, this theory has drawn renewed interests in three 
dimensional gravity in various view points, one of which is along AdS/CFT correspondence. In 
particular, it was shown that the combination of curvature squared terms is consistent with and, 
in fact, determined by the holographic c-theorem, which is a specific incarnation of AdS/CFT 
correspondence [31]~[33- This insight leads to the extensions of NMG to even higher curvature 
theories. Other directions of the exploration of NMG include various (charged) black hole so- 
lutions [2][9]~[l2]) supersymmetric extension [5] [6], generalization by including the gravitational 
Chern-Simons term 00], the computation of correlations and anomaly via AdS/CFT correspon- 
dence p9j, the appropriate extension of Gibbons-Hawking boundary term [39j, the Hamiltonian 
analysis of NMG [55j, etc. 



In this paper, we consider the simplest version of NMG of which action is given 




where r/ and a take 1 or —1, and K \s a, specific combination of scalar curvature square and Ricci 
tensor square defined by 



Our convention is such that m? is always positive but the cosmological constant has no such 
restriction. The equations of motion(EOM) of NMG are given by 




(1) 




(2) 



7? aC^ 



1 



2w? 



= 



(3) 



where 




SR^aK + 2 (4V2i?;.. - V^VuR - g^^u^^R) ■ (4) 



There is an interesting relation between the scalar K and the tensor Kf 



g' 




(5) 



^We have introduced rj for the various sign choice of terms in the action. 
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Another useful form of the action which is equivalent to ([T]) is given by 



S = — 5- / S'x\p^ 



2k 



where fuiy is an auxiliary symmetric tensor field related to 'Shouten' tensor [56] 



and its trace / = g^^ fau- 



2.2 New type black holes 



5, 



(6) 



(7) 



Among various black hole solutions known in NMG: BTZ black holes, Warped AdS^ black holes, 
New type of black holes, Lifshitz type black holes, etc [I] [9] [57], our main interests in this paper are 
the so-called (static) new type of black holes which exist only when a = 1 and l'^ = 1/m?. These 
black hole solutions contain two parameters though the nature of these parameters is obscure at 
the present. One of the motivations for our study is to improve this situation. The metric of new 
type black holes is given by 



ds' 



+ br + c)dt^ + 



dr"^ 



(r^ + 6r + c) 



2m? 



2 



(8) 



with outer and inner horizons at r± 
given by 



i(— 6 lb — 4c). The scalar curvature of this metric is 



2b 



which shows us that there is a curvature singularity at r = when 6 7^ 0. New type black holes 
may be classified according to signs of parameters b and c. This classification may be tabulated by 
introducing a new parameter g = r_/r_|_, as follows: The other sign combinations of b and c do not 
represent black holes. 





b and c 


r+ and r_ 


q = l 


6 < 0, (6^ = 4c ) 


r+ > , r_ > (r_|_ = r_) 


< g < 1 


6 < 0, 


r+ > , r_ > (r+ > r_) 


g = 


6 < 0, c = 


r+ > 0, r_ = 


-1< g < 


6 < 0, c < 


r+ > 0, r_ < (r+ > r_ ) 


q = -l 


6 = 0, c < 


r+ > 0, r_ < (r+ = — r_) 


q<-l 


6 > 0, c < 


r_|_ > 0, r_ < (r+ < r_ ) 



Note that new type black holes for q = —1 are nothing but non-rotating BTZ black holes. 

Now, let us list some quantities of new type black holes. Firstly, the Hawking temperature of new 
type black holes can be read from the surface gravity or the periodicity of Euclideanized action as 

r_i- — r_ \/6^ — 4c , , 
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The mass of new type black holes is identified as a conserved charge by 

M = ^, (10) 

which is also justified by the dynamical approach or AdS/CFT correspondence [M] [35] [58] . Let 
us recall that the Bekenstein-Hawking-Wald entropy of new type black holes is different from area 
law. More explicitly, it is given by 



Ah 

^BHW — 



^'^ at horizon 



where Ah is the horizon area and Ru are Ricci tensors. By computing the entropy on the outer 
horizon, this entropy can be written as the difference between the area of the outer horizon and 
the one of the inner horizon as follows 

Sbhw = ^{A+ - A^) = ^(r+ - r_) = ^V^^c. (11) 

It is interesting to note that the higher curvature effect leads to the inner horizon dependence of 
the black hole entropy. The above quantities satisfy the simple form of the first law of black hole 
thermodynamics as dM = THdSsHW^ which is also consistent with the AdS/CFT correspondence. 



Though some physical quantities are identified and thermodynamic properties are checked, the 
existence of two parameters is still mysterious. It is unclear what is the meaning of another 
combination of parameters different from the mass combination, or what combination at all one 
should take. This problem may be rephrased in another way as follows. 



Usually higher derivative or curvature effects are regarded as introducing additional degrees of 
freedom. This is the reason that the would-be nondynamical three dimensional gravity allows 
propagating massive graviton modes and becomes more interesting. In the form of initial value 
problems of EOM, one needs more initial values to solve the problem with higher derivatives. In 
these viewpoints, the above black hole entropy formula is counterintuitive, since the positive inner 
horizon area means the reduction of the entropy even if it comes from higher curvature effects. More 
distinguishably, there is one parameter family of black hole like solutions, that is, the extremal ones 
(r_|_ = r_), which lead to zero entropy. Therefore, one may envisage the domain of parameters of 
positive inner horizon area as forbidden or at least as disconnected from the case of negative inner 
horizon area. In a later section, we will show some clues on such expectation. 

In order to understand these peculiar properties of the new type black holes, it is desirable to 
investigate these black holes in more detail. In this paper we study scalar perturbations or scalar 
fields on these black hole backgrounds. Specifically, we focus on the QNMs associated with scalar 
fields. It turns out that QNMs give some clues about the nature of parameters. 



3 Exact QNMs of New Type Black Holes 

The EOM of scalar fields on new type black hole backgrounds is written as 

V^^'-m^^' = 0, (12) 
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where m is the mass of the scalar field. To obtain the QNMs of scalar fields on the new type black 
hole background, we recall that the metric of this black hole is given by ([8]). For our convenience, let 
us take L = 1 in the following, which can be recovered by dimensional reasoning. By the separation 
of variables ^' = i?(r)e*'^*^*^'^, one obtains the radial equation of the the above EOM. Through the 
change of variable 

_ r — r_|_ 



(13) 



it becomes 



R"{z) + 



1 



1 



+ 



z z — 1 

1 



+ 



Z- Zo_ 

{z - l){z - Zo) uj'^ (z 



z{z - l){z - Zo) 



ZO 2 

-m 



(z-1) 



Riz) = 0, (14) 



z{zo — 1)2 {z — 1) {z — zq) 

where zq = r^/r^ = 1/^ and ' denotes the derivative with respect to the variable z. 

z°(l — z)^{zq — z)'^%{z) and choosing parameters a, (3 and 7 without loss of 



Defining R{z) 
generality as 



a 



/3 = 1 — \/ 1 + m? and 7 = — V-^o 



(1 - l/^o)r4 

one obtains a differential equation for T-L{z) in the form of 

n'iz) + 



(15) 



■H"{z) + 



2a + 1 _^ 2/3 - 1 _^ 27 + 1 



1 



Z- Zq 



a — 7 



2a7 - 7^ + |(a + /3 + 7)^ - a^jz 



zo(a + /3-2a/3-/32) 
n(z) 



z{z - l){z - Zq) 



0, 



(16) 



which turns out to be just a Heun's differential equation [59j. Let us recall that the Heun's 
differential equation in the standard form is given by 



-H'^z) + 



V 

- + — 

z z ■ 



+ 



Z- Zq 



n'iz) + 



{X^z - r/) 



Z{Z - 1){Z - Zq) 



niz) = 0, 



with the condition e = A + ^ — i^ — 5 + 1. Now, one can see that ()16p is a Heun's differential equation 
by identifying 



1/ = 2a + 1 , 5 = 2/3-1, e = 27 + 1 , 

r] = zo(2a/3 + /3^ - a - /3) + 7^ + 2a7 + a + 7 , 

{A,^} = {2a + /3 + 7, /3 + 7} or {/3 + 7 , 2a + /3 + 7} . 



(17) 
(18) 
(19) 

There exists a local solution of the differential equation (116p near z = which is given by a 
combination of local Heun function //(z)'s. By imposing ingoing boundary condition at the outer 
horizon, one obtains the radial function as 



Riz) = z"(l - zfizo - zyU{z) = Ciz"(l - zf{zQ - zrHizo,7j, A, i^, 5\z) . 



(20) 



Note that as r goes to the horizon (z — t- 0), R{z) is reduced to which represents the ingoing 

mode since a = -n — tt~\ — ■ 
(1-1/20)?-+ 
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Though solutions for the scalar field perturbation on new type black holes can be written as a 
local Heun function, one cannot extract useful information from it since Heun functions do not 
have a closed form in most case and their properties are not well understood, yet. Moreover, 
one should impose the boundary condition at the spatial infinity for the complete specification of 
QNMs. However, in the special cases of new type black holes, we can find QNMs by reducing Heun 
functions to more familiar ones. 



First of all, when q = —1, new type black holes are just the non-rotating BTZ black holes and 
QNMs were already obtained in terms of hyper geometric functions, which can also be derived in 
our approach by the reduction of Heun functions to hypergeometric ones. When q = 0, the radial 
function becomes a (reduced) confluent Heun function. In this case QNMs can also be obtained 
exactly using the connection formula of the (reduced) confluent Heun functions. When q = 1, zero 
modes can be obtained. Now, let us elaborate on these things. 



Firstly, let us consider the q = —1 case, i.e. r_ = — r_|_. The radial function is given by 

R{z) = Ciz^il - zf{-l - zyH{-l, rj, A, C, 2a + 1,2/3- l\z) , (21) 

where rj is given by Eq. ([TS]) while {A, ^} is chosen as {2a + /? + 7, /3 + 7} without loss of generality. 
It has been known that a transformation from the Heun function to the hypergeometric function 
exists when special relations among parameters are satisfied (see Appendix A) [60]. As one can 
check that parameters in the case of zq = —1 satisfy the conditions for the transformation, the 
above solution in terms of Heun function (I2ip . can be written in terms of hypergeometric one. 
Using this transformation, the solution of the radial equation can be written in the form of 



R{z) = C[z''{l-zf{l + z) 



2i^i I a + 



/3 + 7 



a + 



/3-7 



2a + 1 



Az 



(^+1)2 



(22) 



Inserting the values of a, /?, 7 with coordinate change to w = Az/{z + l)'^ , one can see that this is the 
same function already obtained for radial function on BTZ black holes [61]. Therefore, imposing 
the boundary condition (Dirichlet one) at the spatial infinity, QNM frequencies are obtained as 
follows: 

■ ' (23) 



n 



+ -(1 + yjl + m?) 



where n = 0, 1, 2, • • • . 



Secondly, we consider the case (? = 0, i.e. r_ = 0. This case can be obtained by taking a limit 
zq ^ Qo with r_ = rj^/zQ in the Eq. (jl4p . One can observe that this is the limit for the (reduced) 
confluent Heun's differential equation [62]. Therefore, Eq. (|16p is reduced to a (reduced) confluent 
Heun's differential equation 



n"{z) + 



2a + 1 2/3 
+ — 



1 



z z — 1 

with the identification among parameters 



n'{z) 



h^z - rj) 

z{z-l) 



niz) = 0, 



(24) 



a 



lim a 

20 —^00 



UJ 



7 



lim — 

zo->oo Zo 



r.2 ' 



rj = lim — 



7 



a - /3 + 2a/3 + /3^ (25) 
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The local solution of this equation near z = is given by a linear combination of the (reduced) 
confluent Heun function HcizYs. Imposing ingoing boundary condition at the outer horizon, the 
radial function can be represented by 



Riz) = C2z"(l - zfHc 0, 2a, 2/3 - 2, -7^, (/? - 1)2 + f 



(26) 



In order to investigate the behavior at infinity, we can use the connection formula (see Appendix A) 
and obtain relevant part of the radial function for imposing the boundary condition at the infinity 
as 

where ( is determined by 

+ (3 - 2a - 2/3)C + 7^ + - 3a + 2a/? - 3/? + 2 = . 



(28) 



The Dirichlet boundary condition at the spatial infinity {z = 1) requires that the coefficient given 
by gamma functions vanich. This gives the following QNM frequencies: 



^ + {n + \/l + m?){n + ^/l + im? + 1) 
2(n + y/l + m?) + 1 



(29) 



When n is very large, the asymptotic form of QNM frequencies is given by 



(30) 



Finally, let us consider the extremal case (7 = 1, i.e. r_ = r_|_ with denoting the horizon radius as 
ro = r_ = r+. In this case the change of variable in the Eq. ()13p is meaningless. Instead, let us set 
the variable, z = ^-^^ and denote R{z) = z"(l — z)^Q{z). Then, the radial equation in terms of 
Q{z) can be written in the form of 



z 1 — z 

by choosing parameters as 



z{l - z) 



I - Z\(J_^^^ „ 2 q2 
^0 4 



Z 



3 



a^ - - 2a/3 



g(z) = 0,(31) 



a 



2ro 



When u; = 0, it becomes the hypergeometric differential equation. In this case, one obtains the 
solution as 



R{z) 



ciz"(l - zf 2FA- — + a + p, — + a + p,2a + 2 



ro 



+C2Z-''-^{l-zf2Fl 



ro 

— - a + /3-l, — - a + /3-l, -2a 
ro ro 



(32) 



Considering another choice of change of variable and an appropriate ansatz for radial function. 



e.g. z 



r-ro 



and R{z') = z"^ (1 — z'y e^^ G{z') , the radial equation can be reduced to a 
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confluent Heun's differential equatioio. Because of the absence for the connection formula among 
the confluent Heun functions at singularities, it is hard to find exact QNMs for this case. On the 
other hand, there were some other attempts to find QNMs analytically by perturbation method for 
Heun's differential equation [64] [65], 



4 Asymptotic QNMs in New Type Black Holes 



In this section, we present the asymptotic QNM frequencies of scalar fields on new type black 
holes. To obtain these, we use the so-called Stokes line method which matches the approximate 
solutions of the radial equation by the analytic continuation of the radial coordinate to the complex 
plane |49j . It is sufficient to consider the asymptotic solutions near the infinity, the origin and the 
event horizon. After the approximate solutions in each relevant region are obtained, they are 
matched with appropriate boundary conditions. Through these matchings, one can read off the 
QNM frequencies of scalar fields on new type black holes. 



Firstly, we need to solve the scalar perturbation equation (|12p in each region. For this purpose 
it is convenient to introduce the tortoise coordinates x as follows 



dx 
dr 



1 



(r — r+)(r — r_) 

By choosing the integration constant such as x{r = 0) = 0, one obtains 



1 



In- 



r+ 



r — r. 



+ Xq 



where xq is given by 



Xq 



— InT- 



2k^ \r+ 



(33) 



(34) 



(35) 



It is interesting to note that r_|_ should always be taken as a positive number to form black holes 
whereas r_ is allowed to be any value less than or equal to r^. When r_ < 0, there is just one real 
horizon. In the viewpoint of asymptotic QNM computation, r_ = is a special case and should be 
treated separately, whereas the case of r_ < can be done uniformly. In each region, the tortoise 
coordinate x behaves like 

r ^ 
r — )■ r+ . 
r — )• oo 



r_|-r_ 
~ < —CO 

Xq - 



By denoting <I>(x(r)) = y/rR{r), the radial equation of scalar field equation can be written as 



$(x) = 0, 



(36) 



where the potential term is given by 



U{x) 



1- — 1(1 

r 



r 



/i^ + m^r^ + - ^3r^ — (r+ + r_)r — r+r_^ 



(37) 



^It is equivalent to the generalized spheroidal wave equation [63] . 
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The potential term in each region becomes 

' ^+^- (',,2 _ r+r- \ 

U{x) ~ { 



■2 = 



r — > r_|_ (x — >■ — oo) 
r — >■ oo (x — >■ xo) 



In each singular point region, one can obtain the radial function, as follows 



2 2 

C+-^27rw(xo - x) Jjoo (w(xo - x)) 



+ C_-v/27ra;(xn - x)J-j^ (uj(xo - x)) , 

2 



oo {x 



xo) 



(38) 



r r+ (x — )■ — oo) 

where ^^,C± and D± are complex constants and J,i{ujx) and J, joo (i^(xo — x)) represent first 

2 2 

kind Bessel functions. The Dirichlet boundary condition near the spatial infinity implies that this 
radial function should vanish there. Therefore we should take C_ = 0. The condition of £)_ = 
comes from the boundary condition for QNMs near the event horizon ry. $(x) should contain only 
ingoing mode e"^^ there. 

In order to obtain asymptotic QNM frequencies, let us examine the Stokes line which is defined 
by the curve 3m{(jjx) = 0. For the asymptotic AdSo solutions (-D > 4), it has been known that 

the magnitudes of both '2m{ujx) and *He(a;x) are comparable in the asymptotic QNM frequencies, 
which is related to the fact that ux is asymptotically real when x goes to xq, i.e. r — t- oo. Along 
this Stokes line, modes will be oscillatory without any exponentially growing/decaying modes. In 
a neighborhood of the origin, the relation x ~ ^rj^ shows us that 

where p > and n = 0, 1. It describes two lines emanating from the origin, spaced by an angle tt. 
Note that the sign of ux on these lines is alternatively positive and negative. The angle Oq can be 
obtained by 

' for r_ > , 



00 = arg(xo) 



tan 



for r_ < . 



(39) 



In the case of r_ < 0, the branch of the Stokes line emanating from the origin with the angle Oq 
goes to the infinity with positive ux value. The other branch runs to the horizon with negative 
ux value starting from the origin with its angle Oq — tt. The educated guess leads to Stokes line 
depicted in the Figure 1, which results in QNM frequencies consistent with numerical ones given 
in the Appendix B. In r_ > case, the Stokes lines seems to be laid on the real axis. This lines 
do not seem to make sense in matching solutions around singular points. Therefore, we confine 
ourselves to the parameter domain r_ < case only. 
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Figure 1: Stokes line 



Now, we need to match the above solutions along the stokes line on the complex r plane. Using 
the asymptotic expansion of the first kind Bessel function in the region of — t- oo [66J, one 
obtains the expansions of the first and second solutions of Eq. (j38p at some point on the positive 
branch, respectively as 



~ e 



(40) 



4 e 



Matching these two functions gives us a relation 

(joo +J> 



A4. COS 



COXq 



+ cos 



UJXq 



(jc 



(41) 



(42) 



The solution near the origin (|40p is on the positive branch of the Stokes line but the solution 
near the horizon is on the negative one. In order to match the near-horizon solution with the near- 
horizon one, one performs the rotation x' = e~^'^x in the neighborhood of the origin while evading 
the singular origin. By doing this rotation near the origin and then using asymptotic expansion for 
the first solution of Eq. (l38p . the solution ^oix) on the negative branch is represented by 

$o(x) ~ e-'^ {A+e-'^T + A_e''T)e''^'= + e-'^ {A+e-''^ + A_e'^)e-''^'' . (43) 

By matching the above solution with the solution near the event horizon given in the Eq. ()38p one 
can see that ^{x) should be given in the form of e*"^^. Consequently one obtains another relation 

A+e-'— + A.e'— = . (44) 



The existence of nontrivial solutions of two relations ()42p and (I44[) on the coefficients A^ and A. 
implies that their determinant should vanish. As a result, one obtains 



wxo = (n + ^ + ^^vr - ^ In 



/ 771" N 

2cos(:^) 



(45) 
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Note that q is negative since we are dealing with the case r_ < 0. 



Since we know the exact form of xq given in the Eq. (135p . we can obtain the QNM frequencies 
analytically. Substituting k+ = (1 — (?)r+/2 = 2ttTh into Eq. ([35]) . we obtain 



1 



XQ 



AttTh 



Rewriting j in terms of q as 



Hq) 



1 



1 /• 1 1 , 



yj\q\ 2vrrH 

QNM frequencies of scalar perturbation on new type black holes can be written as 



I^QNM 



H 



7r2 + (ln|g|)2 



47r2rff ln|g| 



1 1 



n 



^ i / 

+ + T^Vl + 



27rTf/ In I (7 1 
7r2 + (ln|g|)2 

27r2rjy 



In 



7r2 + (In \q\f V"" + 2 + 2 + + 7r2 + (In \q\f 



In 



(46) 



(47) 



(48) 



In general, the above QNM frequencies show us that the dependence of overtone number n appears 
in both real and imaginary parts. So Jm(w) and £He(a;) take comparably large values with large 
n. This is consistent with the numerical results in large |r_| value region, some of which are 
represented in the appendix B. From the above formula, one can see that the QNM frequencies are 
linear functions of temperature Th ~ \/M for small values, which tells us some information of 
black holes: At least for small \q\, the overtone number dependence indicates that M ~ (r+ — r_)2. 



In the limit of the non rotating BTZ black hole case i.e. q = —1 QNM frequencies become the 
simple form of 

OJQNM = i^ttTh + ^ + ^ \/l + m?^ + 2Th In 

The exact result of QNM frequencies has appeared in [HT], which is given by the Eq. (|23|) in the 
previous section. For the small value of /x and the large value of the overtone number n {\^\ <^ n), 
the real part of the QNM frequencies from the Stokes line method approaches to 2Th In 2 which is 
much smaller than the imaginary part, AttuTh- Therefore we may regard the QNM frequencies of 
the BTZ black hole as a nearly pure imaginary number which is consistent with the asymptotically 
reality condition ujxq for small fi and large n. 



2 cosh 



2Tzj 



(49) 



5 Summary and Conclusion 



In this paper we obtained QNM frequencies of scalar perturbation on static new type black holes 
in NMG. Though scalar field equations were solved exactly in terms of Heun function which is 
the solution of differential equations with four regular singular points, one cannot obtain much 
information from this representation since most properties of Heun functions are yet unknown. 
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Specifically, the so-called connection formulae among singular points are essential for the purpose 
of obtaining QNM frequencies by imposing boundary conditions at different singular points. Never- 
theless, in some special cases, for instance in the BTZ limit, one can reduce Heun functions to more 
tractable ones and then obtain exact QNM frequencies analytically. However, since such reduction 
does not happen in most cases, we have adopted the Stokes line method to obtain asymptotic 
QNM frequencies analytically, which is supported by numerical results for the parameter domain 
— r_ > r+. 

These results reconfirm the mass formula of new type black holes and give some clues about 
their additional parameter. A natural combination of black hole parameters appears in ujqnm is 
the form of Q = ln|g| = ln|r_/r_|_|. Moreover, our results for small Q = ln|q| leads to ojqnm ~ 
Thti ~ n{rj^ — r_), which implies M ~ {rj^ — r_)^, since the overtone number, n, dependence 
of the imaginary part of asymptotic QNM frequencies is usually given by the mass of black hole 
backgrounds, or Choptuik scaling parameter. This is consistent with results from other approaches. 
Though this does not reveal the nature of additional parameter completely, it gives some clues about 
its nature. As alluded in the introduction, we have showed that there seems disconnection between 
r_ < and r_ > during our analysis. That is to say. Stokes line method or numerical one breaks 
down when r_ = 0. More explicitly, two methods are inapplicable in this case because xq cannot 
be defined and the potential of the radial equation in Eddington coordinates is unbounded. This 
may not be just inadequacy of methods to the problem as was argued by the entropy reduction in 
the r_ > domain. 

There are still many issues to be addressed in the future. First of all, it is very desirable to 
study AdS/CFT correspondence for QNM frequencies in new type black holes. Since new type 
black holes have different metric fall-off tails at the spatial infinity from BTZ black holes, the dual 
CFT may be different. Nevertheless, it is still perplexing that there is an additional parameter 
in the dual CFT picture. Thermal two point functions of perturbation operators in CFT can be 
determined more or less uniquely only by conformal symmetry and depends only on the conformal 
weights, which have been identified as black hole mass or temperature. It is unclear how to identify 
the additional parameter in the dual CFT. Maybe new type black holes are not good quantum 
objects allowed in quantum gravity. 

Secondly, though we presented exact QNM frequencies for r_ = and r_ = — r_|_, these cases 
should be studied more carefully, since the potential in tortoise coordinates can have negative values 
when < — r_ < r_|_ or r_ > 0. This can be referred to as the necessity of further study about 
stability of new type black holes. During our analysis, we have assumed that new type black holes 
are stable under perturbations. However, this is unverified assumption, though it is very plausible 
that new type black holes are stable near the parameter region corresponding to BTZ black holes 
which are stable since those can be embedded as Bogomol'nyi-Prasad-Sommerfield (BPS) states in 
the supersymmetric extension of NMG. 

Thirdly, it is worthwhile to consider the problem on the quantization of the black hole entropy with 
the asymptotic QNMs. There have been many investigations concerned with this |52j |53j [67]^|72j. 
According to the work of |71j|72j. we can calculate the entropy spectra of new type black holes 
from asymptotic QNMs (i.e, for large overtone number). When q = —1 and q = 0, the entropy 
spectrum seems to be equally spaced. However, when q < 0{q ^ —1), the spacing of the entropy 
spectrum depends on the q factor, even though the entropy spectrum is still equally spaced for fixed 
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q value. The different aspect of the entropy spectrum may be caused by the different behavior of 
the asymptotic QNMs with both highly oscillating and damping modes. Some similar phenomena 
are observed for large Schwarzschild-AdS black holes in D > 4. Further investigations are needed 
to understand the entropy spectrum using the asymptotic QNMs. 

Finally, though numerical results support Stokes line or monodromy method in general, the 
rigorous mathematical justification of this method is still lacking. More concretely, it is unclear 
how to estimate errors in Stokes line method and how to systematically compute the next orders 
in this approach. Usual convergent or asymptotic scries approach to QNM frequencies has definite 
error estimates and can be investigated systematically, which leads to numerical computation of 
QNM frequencies. Therefore, it is also very interesting to investigate Stokes line method rigorously 
and verify its validity in general. 
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A The reduction of the Heun function to the Hypergeometric 
function 



In this appendix some mathematical identities are presented, which are useful in the main text. 



A.l Heun to hypergeometric 



The reduction formulae from the Heun function to the Hypergeometric one by polynomial transfor- 
mations are given in mathematical literatures, for instance see [60] . In summary, when four singular 
points, (0, l,zo,oo) of the Heun function satisfy harmonic or equianharmonic conditions there are 
polynomial transformation from Heun function to Hypergeometric one. The harmonic condition 
means the three points (0, 1,20) are collinear and equally spaced, which is relevant to our case. 
In [60j . zq = 2 is taken as the canonical value with the polynomial transformation R{t) = t{2 — t), 
see the Eq. (3. 5. a). The reduction formula in the case of zq = — 1 can be obtained by choosing a 
Mobius transformation from z to t as 

t 



2z 



z + 1 

(^^^\y2 and (3. 5. a) actually means two equations, one can see that 



Noting that R{t) = t{2 - t) 
there are two possible ways to reduce Heun functions to Hypergeometric ones as follows. That is 



H{-1, r?. A, i,u,5\z) = [z + 1)- 2^1 ( ^, ^ ^ ^ 7 



4z 



[z + lf 



(50) 



with relations among parameters as 



(51) 



or the similar formula with A and ^ interchanged. This formula explains the reduction of radial 
equation from new type black holes to BTZ ones when r_ = — r-)_. Here, our convention of the 
Hypergeometric function is 







l + a + h 



w 



1 



c\ d 
J dw 



ab 



■w{l — w) 



2Fi{a,b,c\w) . 



(52) 



A. 2 Confluent Heun to hypergeometric 



The (reduced) confluent Heun's differential equation is given by 

The general solution of this equation near z = is given by a linear combination of the (reduced) 
confluent Heun function Hc{zys 

n{z) = ciHciO, b, c, d, e\z) + C2z'''Hc{0, -b, c, d, e|z) (54) 
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where e = i{l - (6 + l)(c + 1) - 2e}. 



There is a connection formula which is given by 

Hc{0,b,c,d,e\z) = Cl Hc{0, c, b, -d,e + d\l- z) 

+^^ r(l + c|a)r(b-c7) ^^ - ^r^HciO, -c, b, -d,e + d\l-z), (55) 

where 

(^■^ + ^l-b-c)C-e-b-c+^ = 0, 

(7^ + (1 - 6 + c)(7 - e - 6(c + 1) + ^ = . 

The above formula means that the confluent Heun function about z = can be connected with 
some combination of the two solutions about z = 1 by analytic continuation. 

When d = 0, the reduction formula from the (reduced) confluent Heun function to the hyperge- 
ometric one is given as follows: 

HciO,b,c,0,e\z) = 2Fi{A,B,C\z) , (56) 

where 



{A,B} ^ 



6 + c ± + c2 - 4e + 1 + 1 6 + c T Vfe^ + c2 - 4e + 1 + 1 



C = 6+1 



A. 3 Hypergeometric function and connection formula 

The standard form of the hypergeometric differential equation is given by 

The solutions near z = is given by 

F{z) = (a, b, c\z) + z^~''2Fi{a - c + l,b - c + 1,2 - c\z) . (58) 

The connection formula from the hypergeometric function, 2Fi{a, b,c\z), about z = to another 
one about z = 1 is given by 

2F,{a,b,c\z) = {l-zr—'^^^^^^^^^-^^2Fi{c-a,c-b,l-a-b + c\l-z) 

r(a)r(6) 

+ ^&t^-±±^^F,{a, b,l + a + b-c\l-z). (59) 
r(c — a)T[c — b) 
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B Numerical results 



In this appendix, we present some details about the numerical computation of QNM frequencies 
following |48j [74j [75] . New type black holes with ingoing Eddington coordinates (f; = t + x) are 
given by 



f{r)dv'^ + 2dvdr + r 



2 Jj,2 



/(r) = (r — r+)(r — r_) . 



(60) 



Scalar fields on these black hole backgrounds may be decomposed as ip = ;^$(r) g^'^^e*^'^ under the 
convention taken in this paper. Note that in this convention ujqnm should have positive imaginary 
part. The above separation of variables leads to the radial equation of massless scalars as 



/(r)$"(r) + (f'{r) + 2iuj) $'(r) - V{r)^{r) = , 



d 
dr 



(61) 



where 



V{r) 



r_|- + r_ 
4r 



+ 



fi — r^r^ 

47-2 



By introducing new 


variable y = 


1/r, 














where y+ = l/r_|_, y_ 


= l/r_ and 








(y 




y+y- 




t{y) 


2 


y{y 




^ y+y- 




u{y) 


= (y-y^ 


OF(r 


Series solution ansatz 


is taken as 





i(y) d 



y-y+dy (y - y+) 



r$ = 0, 



y+ + y- 2 , o- 2 



Hy) = ^a^(y-y+)^ 



^=0 



which gives us recursion relation of a£ in terms of the coefficients of mode expansion s{y),t{y),u{y) 



ae 
Pi 



A:=0 



l)so + Ito = 2y\l{^K+ + iuj) 



These formulae can be implemented easily in mathematica [76] . which leads to the following results 
in the case of — r_ > r i as follows. 
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r+ 


r_ 


fJ- 






r+ 


r_ 






1 


-10 





18.0486 20.9188 i 




100 


-200 





108.272 + 284.046 i 


1 


-10 


1 


18.5137-h 20.5174 i 




100 


-200 


1 


108.282-h 284.044 i 


1 


-10 


2 


19.5422+ 19.6755 i 




100 


-200 


2 


108.312-F 284.036 i 


1 


-10 


3 


20.7523-h 18.7317 i 




100 


-200 


3 


108.363 + 284.024 i 


1 


-10 


4 


22.0311-F 17.7672 i 




100 


-200 


4 


108.433 + 284.006 i 


1 


-10 


5 


23.3505-1- 16.7959 i 




100 


-200 


5 


108.524 + 283.984 i 




r+ 


r_ 








r+ 


r_ 






10 


-20 





10.8272 + 28.4046 i 




10 


-10 





0.00000-F 20.0044 i 


10 


-20 


1 


10.9274 + 28.3797 i 




10 


-10 


1 


1.00022 + 20.0022 i 


10 


-20 


2 


11.2203 + 28.3069 i 




10 


-10 


2 


2.00045-F 20.0022 i 


10 


-10 


3 


11.6854 + 28.1916 i 




10 


-10 


3 


3.00068 + 20.0022 i 


10 


-10 


4 


12.2947 + 28.0413 i 




10 


-10 


4 


4.00091 + 20.0022 i 


10 


-10 


5 


13.0192 + 27.8644 i 




10 


-10 


5 


5.00113 + 20.0022 i 



One can see that there are very smah variations of the numerical results when /x values are changed, 
in particular for large r_|_ and r_. Of course, this numerical results are just taken for the low lying 
modes and the number of interaction are taken up to 50 times i.e. £ = 50, but it strongly indicate 
the highly damped mode obtained by Stokes line methods is reliable at least — r_ > r+. 
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